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Abstract 

The question of whether the zero viscosity limit v — > is identical to the no viscosity 
v = case is investigated in a simple shell (GOY) model with only three shells. We 
find that it is possible to express two velocities in terms of Bessel functions. The third 
velocity function acts as a background. The relevant Bessel functions are infinitely 
oscillating as v — > and are not analytic as functions of v at the point v = 0. We also 
mention a perturbative method which may be used to improve the model. 
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1 Introduction 



In the Navier-Stokes equation the viscosity v plays an important role. It is an interesting 
question of principle whether the limit of vanishing viscosity is given by the solution to 
the same equation with v identically zero. In other words, is the limit v — > smooth, or 
is it non-analytic? An answer to this question may be of interest in various branches of 
physics, for example in cosmology, where the history due to the expansion of the Universe 
may contain different viscosities and the behavior during transitions from one v to another 
may not be analytic. 

In a numerical study of decaying turbulence some evidence was found in [1] that the 
small viscosity limit is highly non-trivial and does not conform to the naive expectations. 
Of course, numerical data do not necessarily allow an extrapolation simulating the limit 
v — > 0. Therefore it would be desirable to have some explicit mathematical expression for 
the velocity field which allow this limit to be investigated. 

It is well known that a direct mathematical study the Navier-Stokes equations for high 
Reynolds numbers is not an easy matter. Therefore there is some motivation for looking for 
a model of the hydrodynamics equations where an analytic approach may be more hopeful. 
To this end one may think of the shell (GOY) model of Gledzer, Ohkitani and Yamada[2]. 
Many properties of turbulence, especially those related to energy transfer and the small 
intermittency effects, have been understood from the numerical studies of the shell model. 
For a review of the applications, see ref . [3] . 

The model is formulated in terms of Fourier space velocity variables u n (t), and the 
dynamical equations are given by 



Here k n = r n , and / is an external force acting on shell number no and n is less than some 
maximum number N. Usually this equation is studied numerically with a large number of 
shells. Also, the usual scaling law of Kolmogorov k^ 1 ^ 3 appears when the maximum number 
N of shells go to infinity [3]. 

There is not much hope that the shell equation (1) can be integrated in terms of standard 
mathematical functions. Therefore the question of whether the solution of (1) in the limit 
^ — > is identical to the solution of (1) with v = would still be subject to numerical 
extrapolation. 

However, the probability of obtaining explicit solutions for the u' n s may increase if the 
number of shells is small. In this paper we shall show that a simple truncated three-shell 
model actually allows a solution in terms of Bessel functions, provided the parameter 5 takes 
a somewhat special value. These Bessel functions turn out to be infinitely oscillating in the 
limit v — > 0, so that as v = is approached, they can take many different values, and in 
general the v — > dynamics does not correspond to the v — case. In some cases it is 
possible to "renormalize" the bad behavior so as to obtain a smooth limit for small v, at the 
cost of trading a quantity which is infinite at v — with an initial velocity. 

In section 2 we introduce the three-shell model. As a preliminary we solve it for v = 0. 
Then in section 3 we include viscosity and solve the model in terms of Bessel functions. In 
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section 4 the limit v — > is investigated. Section 5 contains a discussion of the large time 
limit, and section 6 discusses the relatively simple case where the shell distance is very large. 
The Bessel functions are then approximately replaced by trigonometric functions, which 
are still oscillating infinitely near v — 0. In section 7 we discuss a possible perturbative 
improvement of the model, and we conclude in section 8. 



2 The three-shell model 

We shall now consider a simple GOY-model with only three shells, and with the forcing term 
acting on the first shell. Furthermore, we consider the case where 

5 = 1. (2) 

This value means that the dynamics is intermediate between two dimensions, where 5 = 
1 + 1/r 2 , and three dimensions, where 5 = 1 — 1/r. The choice of 5 is connected to invariants 
like enstrophy or helicity invariants, as was first pointed out by Kadanoff et al. [4] . See also 
the review [3]. The general equations (1) now reduce to the following three equations 

(it + Uk *) = - lk ^ t >^ t ) + /' ( 3 ) 

and 

( j t + vk\ J u\{t) = ik lU3 (t) Ul (t), (4) 



as well as 



j f + vk£\ ul(t) = 0. (5) 



The last equation simplifies because of S = 1, since otherwise the right hand side would 
contain the term ik\{l — 5)uiU 2 - 

For future reference we first solve these equations with zero viscosity, v — 0. We impose 
the boundary conditions 

ui(0) = and u 2 (0) ^ 0. (6) 

Eq. (5) can be trivially solved, 

u 3 (t) = C 3 , (7) 

where C 3 is a (complex) constant. Because of this, the velocity field u 3 (t) will act as a 
background field for the other velocities U\ and u 2 . Using this, eq. (3), which now reads 

d -^f> = -ih C 3 u 2 (t) + /, (8) 

can then be integrated, 

u *(t) = ft- ih C 3 [ dt' u 2 {t'). (9) 
J o 
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Inserting this solution for U\ in eq. (4) with u 3 replaced by C 3 we obtain 

+k 2 1 \c 3 \ 2 v*(t) = ik 1 c 3 r- (io) 



^M*) , ,2,., ,2 . , 



dt 2 

This second order differential equation has a simple oscillating solution which can be inserted 
in eq. (8) in order to obtain u±. Using the boundary conditions (6) we simply get 

2iC iC f* 

u *(t) = — f a sinjfei|C 3 |t, u%(t) = 2a cos h\C 3 \t + , and u* 3 (t) = C 3 . (II) 

Ml " h\C 3 \ 2 

Here ct is an integration constant which is related to u 2 at t = 0, 

Um = 2a+ k^- (12) 

It should be noticed that although the forcing was coupled to the equation (3) for u 1: in the 
solutions (11) the constant / occurs also on the second shell. 

3 The three-shell model and viscosity 

We now consider the basic equations (3)-(5) with non-vanishing viscosity, v ^ 0. Again, eq. 
(5) can be solved trivially, 

u 3 (t) = C 3 e^ fc a 2 *. (13) 

Again this field will act as a background field for the two other velocities. Inserting eq. (13) 
in (3) we can integrate to obtain, 



u*(t) = e""*?* / dt' e^'fti') - ih C 3 e' uk " 1 [ dt' e -^I" fc ?)*' u 2 (t'). (14) 

JO JO 

Inserting this in eq. (4) gives 
(| + «SW 

= -kl\C 3 \ 2 e-"«+*»" /' A'e-""-? -'?>•' (■*(«') + iC 3 kj e -«+*»>' f dt'e""' 1 ' f'(t'),(lf,) 

JO JO 

and by differentiation we obtain the following second order differential equation, 



^ + v{k\ + kl + kl)^ + {v 2 k\{k\ + k 2 3 ) + kl\C 3 \ 2 e~ 2 ^)u* = ihC^e-^. (16) 

It should be noticed that formally eq. (16) reduces to (10) by taking v — 0. The main point 
of this paper is that this is not true for the solution of (16), unless time is very small. 
In order to solve eq. (16) it is convenient to introduce the function S(t), 

S(t) = e ^l+kl+ki)t/2 u ^ t y (17) 
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Eq. (16) then gives 

^ + + kl + klf + Skl{k\ + A; 2 ) + fc 2 |C 3 |V 2 ^ S 

= ihCztity^ 2 ^-^ 1 (is) 

The homogeneous equation corresponding to / = is of the Bessel type, with the solution 
where 



, 2 v '(fc 2 + *£ + A; 3 2 ) 2 - 4A; 2 (A; 2 + kl) = (kl - k\ + kf). (20) 

Since k n = r n we see that for large r the index a approaches 1/2. 

The inhomogeneous equation (18) can be solved by standard methods from the knowledge 
of the solution (19) of the homogeneous equation, 

S (t) = ttt^-. ikiCs (j-a(z) fdt'J a {z')g{t')f{t') + J a {z) f dt' J_ a {z')g{t')f{t')\ , 

zuk^simra \ Jc Jt J 

(21) 

where c and d are arbitrary constants, and where 



g( t ') = e -Hkl~kl~ki)t>/2 (22) 

and 



and z' is the same as z except that t is replaced by t' . 

It should be noticed that the argument of the Bessel function (23) is not nicely behaved 
as a function of v. Therefore the formal coincidence of eqs. (10) and (16) for v = is not 
reflected in a simple manner in the Bessel solution (19). 



4 The limit v -> 

In view of the last remarks we shall now study the limit v — > 0. To simplify matters, we 
start by looking at the homogeneous solution 

u\(t) = AJ a (z) + BJ_ a (z), (24) 

where A and B are integration constants, and where z was defined in eq. (23). If we maintain 
the boundary condition -Ui(O) = 0, it follows from eq. (4) that 

= 0. (25) 

t=o 



^ + vkluX 
dt 1 1 
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From the asymptotic form of the Bessel function valid for v — > 0, i.e. z — > oo, 

J ±a (z) « J— cos(z T an/2 - tt/4), (26) 

V 7IZ 



and from the explicit expression for z in (23) it follows that behaves like s/v times a 
cosine. The time derivative of also behaves like ^/v times a sine, so therefore vk\ix\ 
(which behaves like z/ 3 / 2 ) is subdominant relative to the time derivative of u-z when v — > 0. 
Therefore we can replace the boundary condition (25) by duijdt = 0. This allows us to find 
the following relation between the constants A and B in eq. (24), 

B = -A 2 4 -{, (27) 

where we used the asymptotic form (26) for the Bessel functions. Introducing instead of 
A the initial value 112(0), which is assumed to be independent of u, we now get by use of 
addition theorems for trigonometric functions 



«$(o) 



nh\C 3 \ -(^ + af-jO (2g) 



2 sin(7ra) 



This result allows us to trade the integration constant with the initial velocity, and by use 
of (26) we have 

u*(t) » u*(0)e'^ +k ^ k ^ 2 cos (l - e'^fj . (29) 

Again we used addition theorems for trigonometric functions. This result agrees completely 
with the result (11) for v — with / = if the expansion exp(— vk 2 i) ~ 1— vk\t is performed. 

In this case it is thus possible to absorb the bad behavior for v — > in the "renormaliza- 
tion" (28) of the constant A in terms of 112(0). The price to pay for this is that the badly 
behaved constant A (for v — > A blows up with infinite oscillations) is traded with an 
initial velocity assumed to behave nicely as a function of v. Therefore, if one does numerical 
simulations it would look as if the limit where v is decreased more and more does not lead 
to definite convergent results, unless the "renormalization" discussed above is performed at 
each step of the limiting procedure. 

The basic reason for the violent behavior of A is that the Bessel function (19) is not 
analytic as a function of v in the point v — 0. Here the most singular behavior is 

%^ « -mr* - ftf ^' * -5 - 3 + (£) ■ <3o) 

This shows that the solution of the differential equation (16) does not behave in a simple 
way for small v. 

It is not always possible to "renormalize" the bad behavior away. The function 

u\(t) = ^(0)e-^ +fe 2 2+fe 3 2 )*/ 2 J a (^|C 3 |/(^3 2 ) e-^/Uhld/^kl)) (31) 
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is a solution which rapidly oscillating for v — > 0, and no "renormalization" trick can remove 
this behavior. We thus see that the simple three-shell model can produce a variety of results 
which look rather non-trivial. 

Instead of imposing some initial velocity for the u\ and u 2 fields, we can proceed in a 
more physical way by assuming that initially u\ = u 2 = up to the time to- Then some 
external agency applies a force for a limited time from to to ti, so after the time t± the 
solution of the homogeneous equations of motion emerges, as can be seen from eq. (21) with 
f(t) inside the integrals. However, if we write this solution as 

S(t) — A J a (z) + B J-a(z) for t > ti, (32) 

then A and B are not arbitrary coefficients to be fixed by some initial velocity of Ui and/or 
u 2 . On the contrary, these constants are fixed dynamically by the force and by the initial 
velocity C 3 of u 3 , 

A = - ~!f lCz r dt J- a (z) g(t) f(t)B = 3 f 1 dtJ a {z) g(t) p(t). (33) 

2 v sm7ra Jt 2 v k% sm7ra J to 

Here there is no way of "renormalizing" the constants. If the force acts in a short time, the 
situation is similar to what was discussed above. For example, if we take a delta-function 
pulse, 

f(t) = fo6(t-U), (34) 
we obtain for v so small that git) 1 

s (t) « Z r^ sin ( ^nr ( e - vklt ~ e ' uqu )) W ~ **)■ (35) 

IG3I \ vk 3 v ' ) 

In arriving at this result we used standard addition formulas for the trigonometric functions. 
For v small and time not too large this again gives the v — result, which in the present 
case is arrived at smoothly. 

The situation changes completely if we let the force act for a longer time. Then the 
integrals in eq. (33) in general involve times where an expansion like e^ vk -'- t wl - vk\t is not 
valid inside the integrals. Therefore a result analogous to (35) with a smooth v — > limit 
does not appear if the force is allowed to act long enough. As a matter of fact, the integrals 
determining A and B in (33) are strongly oscillating, as one can see in numerical examples. 



5 The large time limit 

In the case where t — > 00 in such a way that \ jv e~ vk '& — > only the lowest order terms in 
the Bessel functions should be kept, J a oc z a (l + 0(z 2 )), so for / = we have 

u* 2 (t) » L_ e^- 1 ' 2 { 1 - k 3^L e - 2 » k2 A + L + e~^ 2 ( 1 - e'^A , (36) 

2V ; V 4r ( 2 -°) / V 4r ( 2 + a ) I 

where the L's are constants and 

b± = k 2 + k 2 2 + kl± 2k 2 3 a, (37) 
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which means 

b + /2 = kl + kj and 6_/2 = k\. (38) 

The function u\(t) can be obtained most simply directly from eq. (4). The result is again 
expressed in terms of Bessel functions. We can find the asymptotic behavior corresponding 
to eq. (36) by inserting eq. (36) in eq. (4). From (38) we have for t large 



u* 2 (t) » L_ e-^ ( 1 - k 3C*L e -*>%t\ + L e -u{kl + kl)t ( j _ kj\C 3 \ 2 ^ 2vklt | 

4r(2-a) ; V 4r(2 + a) 



The leading term contains e - ^ 2 * as one would expect. It is important that this term is 
annihilated by the operator d/dt + vk\ in eq. (4). By use of eq. (4) the result for t large is 
thus 

, L + kl\C 3 \ 2 (3kj- kl + kp {k 2 +2k 2 )t 

+ 4r(2 + a) j 1 } 

As one would expect, the leading term is simply e~ uk ^- t , which is annihilated by the operator 
d/dt + vk\. This u\ can be reinserted in eq. (3) as a check of the self-consistency of eqs. (3) 
and (4). 



6 Simplified results for r ^> 1 

We shall now mention that in the limit where the shell distance r is very large, our results 
simplify considerably. From eq. (20) we obtain a ~ 1/2, since k\ and k\ can be ignored 
relative to k\, so the solution can now be expressed in terms of trigonometric functions, 



V Trfci I C 3 1 



a ■ I ki\C 3 \ _,,i,2 i -\ I ki C3 _ v u2 f 

A sin ' ' e vfe 3* + 5 cos 1 



z/fcg / y vk\ 

+ inhomogeneous / term. (41) 

From eq. (4) we can then find the corresponding field Ui(t), 



.-(A'r - fc 2 2 ) + fcilCalS e-^*) sin (^§1 

2 V Z/rCq 



£(jfe2 - A; 2 ) - A^ICale-^ 2 * + vk {B) cos f ^r^ 1 ' 



'3 



+ inhomogeneous / term. (42) 

Again we can fix the integration constants by suitable boundary conditions. Also, we see 
that the limit v — * is not well defined due to the non-analytic behavor of the arguments 
of the sine and cosine, which oscillate violently as v is decreased. 



8 



7 Perturbations in 1 — S 



To improve the approach presented above one could try a perturbative expansion in (1 — 5). 



Denoting the 5 = 1 functions found above by uf^ and u 2 u> eq. (5) would change into 



,(°) 



d 



j t + z/A; 2 I u* 3 (t) = - 5) u^u { 2 0) 



(43) 



This is an equation for u% with a time- dependent "forcing" term. It can be solved for u% in 



terms of the unperturbed functions and u y 2 ' . This perturbed u$ should then be inserted 
in eqs. (3) and (4) to give the perturbed u-y and u 2 . Of course, there is no guarantee that 
such an expansion is convergent. 
From eq. (43) we easily obtain 



(o) 



|u 3 | 2 « e - 2 "*3 a * 



If we use eq. (3) to express u 2 in terms of ui, so we have 



(44) 



c s «i o) (t)t4 0) (*)-c?(«r(or« ; (*)) 



(0), 



Mi 



■ vk 2 t 



d\u\ 



(0)|2 



+2^i4 o) i 2 -/< ; -/x 



(0) 



dt 



,(0)* 



This expression can be inserted ie eq. (44) and after a partial integration we obtain 



F3 



-2vk%t 



where we defined 



and 



|C 3 | 2 -(1-5) (e 2 ^*|4 0) (t)| 2 - |mS 0) (0)| 2 ) + (1 - 5)I(t) + (1 - 5)J(t) 

I(t) = 2v(kl-kj) f 
Jo 



dt' e 2 ^' 



(0 



j(0 



(MV) + r(4V)) 



(45) 

(46) 
(47) 

(48) 



It is interesting that the third term in the square bracket in eq. (46) containing the integral 
/ has a definite sign depending on 1 — 5: If 1 — 5 > the sign is positive since ks > k 2 - This 
5 corresponds to helicity (= XX~ l) n k n \u n \ 2 ) conservation for v — f — 0. So through the 
/—term the two first shells give a positive contribution to the energy of the third shell in 
the lowest order perturbation theory, which means transfer of energy from lower to higher 
k's. This is precisely what would be expected in three dimensions. 

On the other hand, if 1 — 5 < 0, the /—contribution is negative, and the energy in the 
third shell is decreased by this effect, as expected in two dimensions where the enstrophy 
(— 2 kn\ u n\ 2 ) is conserved for v — f — 0. 

Of course, in the full expression (46) there are two other terms proportional to 1 — 5. For 
/ = it should be possible to investigate all the terms in (46) numerically by inserting a 
Bessel function constructed from the equation of motion (4) and the solution for u 2 , thereby 
giving u\. It would be interesting to see if the overall sign is +(1 — 5) corresponding to the 
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expected inverse cascade (transfer of energy from shorter to larger scales, i.e. from larger 
to smaller fc's) in "two dimensional" systems (5 = l + l/r 2 >l) and a forward cascade in 
"three dimensions" (5 = 1 — 1/r < 1) where the energy is transported from smaller to larger 
k's. 

It should be noted that if r is large, 5 is in both cases close to 1. Consequently, if the 
three-shell model should attempt to be somewhat similar to the GOY model with many 
shells, this will work best for large separations between the shells. 

8 Conclusions 

In the simple three-shell model we have found the velocity functions in terms of Bessel 
functions with an argument which is not analytic in the viscosity v. In some cases it is 
possible to hide this singular behavior by a suitable "renormalization" . However, this is not 
true in all cases. So the model may never approach the similar model with no viscosity, 
v = 0, even if v — > 

Of course, the reason for the integrability of the three-shell model is that the third shell's 
velocity becomes a fixed background field for u\ and u 2 - Thereby the complexity due to 
the basic non-linearity of the shell model has disappeared or, at best, become rudimentary. 
However, this does not make the model completely trivial, since there is a non-trivial coupling 
between u\ and u 2 . So different Fourier modes do couple, and the coupling makes transfer 
of energy between these modes possible. 

We also discussed a perturbation approach with an expansion in 1 — 5 around 5=1 where 
the three-shell model was originally defined. Although the resulting perturbative change of 
the energy of the third mode is relatively complicated, we identified one term which has 
the expected cascade properties between three and four dimensions. We hope to be able to 
perform an precise analysis of the sign of all the terms later. 

In conclusion one can say that although the three-shell model is of course an immense 
simplification of the GOY model, nevertheless it has features which indicate a non-trivial 
dynamics. 

Acknowledgement : 

I thank Mogens H0gh Jensen for many interesting discussions over several years. 

References 

[1] T. A. Yousef, N. E. L. Haugen, A. Brandenburg, Phys. Rev. E69 (2004) 056303. 
arXiv : astro-ph/0312505 

[2] E. B. Gledzer, Sov. Phys. Dokl. 18 (1973) 216; K. Ohkitani and M. Yamada, Phys. Rev. 
Lett. 60 (1988) 983. 

[3] T. Bohr, M. H. Jensen, G. Paladin and A. Vulpiani, Dynamical Systems Approach to 
Turbulence, Cambridge University Press (1998). 

[4] L. P. Kadanoff, D. Lohse, J. Wang and R. Benzi, Phys. Fluids 7 (1995) 617. 



10 



